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Abstract 

Various implications of new, non-perturbative pomeron inspired 
enhancement of small- x neutrino-nucleon structure functions for 
high-energy neutrino astrophysics are discussed. At x > 1CP 5 
these functions are given by perturbative QCD, while at lower x 
they are determined by a specific generalization of F 2 ep (x,Q 2 ) de- 
scription, proposed by A. Donnachie and P. V. Landshoff (their 
two-component model comprises hard and soft pomerons), to vN- 
scattering case. 

We found that i) such enhancement causes the most rapid growth 
of i/iV-cross-sections at high energies, ii) that pomeron effects may 
be perceptible in the rates of neutrino induced events in future gi- 
ant detectors and iii) that the rate of high-energy neutrino flux 
evolution (due to absorption (CC + NC) and regeneration (NC)) 
on its pass through a large column depth of matter may be sub- 
jected to additional influence of hard pomeron. Solving transport 
equations for the initially power-law decreasing z^-spectra, we have 
evaluated shadow factors for several column depths and spectrum 
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indices. The results are compared with analogous calculations, per- 
formed within a trivial small- x extrapolation of structure functions. 
Hard pomeron enhanced high-energy shadow factors are found to 
be many orders of magnitude lower than those obtained within or- 
dinary perturbative QCD. 

1 Introduction 

As a rule, accelerator experiments provide data of very high quality. But 
such measurements seem to be hardly possible at energies, E > 1 x 10 19 eV, 
even in the very distant future. 

On the other hand, non- accelerator, observational physics gives a good 
chance to get data at ultrahigh energies M. But the accuracy of data 
in this case is worse due to scarcer statistics and lower resolution of de- 
tectors. It is remarkable, that non-accelerator High-Energy (HE) physics 
combines into one the 'telescope' and 'microscope' physics. Really, the 
detected UltraHigh- Energy (UHE) particles are produced in perhaps most 
distant, extragalactic cosmic sources. Hence, it deals with extremely large, 
cosmological, distances. On the other hand, interactions of such particles 
with matter occur at least distances. 

Realization of UHE astrophysics is a complex problem. It poses many 
hard questions, such as: 

1. what types of ultrahigh-energy particles may be observed? 

2. where and how ultrahigh-energy particles are produced? 

3. what detectors are to be used? 

4. what is the expected rate of events in a detector? 

5. does this rate exceed the background? 

Let us try to answer these questions in brief. 
1.1 Particles 

Through the centuries people used only star light, viz. optic photons, for 
the sky observation. Today actually the whole range of e.m. radiation, up 
to TeV photons 0, is involved in measurements. 

Another well-known example of high-energy cosmic radiation is pro- 
vided by Cosmic Rays (CRs). These are mostly protons (if), He, C, N 
and up to Fe ions, which continuously bombard the atmosphere. Energy 
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spectrum of CRs extends from 1 GeV up to perhaps 1 x 10 21 eV; the high- 
est energy event with E ~ 3 x 10 20 eV has been detected by the AG AS A 
array [[|. The nature of incident particles in such events as well as where 
and how they are produced is unknown. Moreover, the average path length 
of E > 6 x 10 19 eV protons in the cosmic microwave background radiation 
(with temperature T = 2.73K) is essentially limited by p + 7 — > 7r + + n 
interactions (GZK cut-off), so that their presence in CRs is a problem. 

Besides traditional photons and CRs, gravitation waves and neutrinos 
are also regarded as possible instruments for HE astrophysics. In this paper 
we concentrate on neutrinos, in particular, on UHE neutrinos with E v > 
1 x 10~ 15 eV. Being neutral, they do not deflect in magnetic fields, hence 
their arrival directions shoot back to the production site. Small cross- 
sections allow these neutrinos to travel extremely large distances without 
absorption. But their production in cosmic sources is difficult (i.e. sources 
are rare) and these sources are distant. It makes the expected flux to be 
small. Moreover, detection of such neutrinos is complicated by the same 
small cross-sections. To compensate for these negative factors, one should 
construct gigantic detectors. 

1.2 Production 

There are many theoretical models predicting large luminosity of cosmic 
sources in high-energy particles, including neutrinos (see e.g. [I]]). For 
example, such fluxes may be produced in supernova explosions, in active 
galactic nuclei, during gamma-ray bursts etc. Energy outburst in these 
sources, including in the form of HE particles, may be really high, but the 
expected flux at the Earth remains low due to large distance. 

Basically, there are two models of HE particle production. First as- 
sumes proton acceleration up to very high energies on cosmic source shocks 
(usually gas accretion onto a massive black hole, presumably residing in 
center of an active galaxy, is suggested), sometimes in relativistic jets point- 
ing to the Earth. Then pions and kaons appear as a result of pp- and/or 
jry-collisions in a source. These mesons decay, giving birth to HE neutri- 
nos. The described scenario is called down-top mechanism. Unfortunately, 
it is difficult to obtain proton energies higher than E ~ 10 17 eV in this 
mechanism; one needs both large shock radius and high magnetic field. 

Another possibility is proposed by the top-down mechanism (for review 
and discussion of upper limits to the expected ^-fluxes see e.g. Ref. || and 
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references therein). It is based on suggestion that some yet unknown, 
extremely massive (with mx ~ 10 14 -j- 10 16 GeV) long-living gauge X- 
bosons decay to the ordinary particles. These bosons allegedly originate 
from topological defects, which in their turn have been produced during a 
hypothetical phase transition in the early Universe. A big fraction of X- 
boson's energy goes to HE particles (neutrinos). These attractive models 
involve new physics beyond Standard Model and seem a bit speculative. 
An experimental verification is needed. 

1.3 Detectors 

HE neutrino detectors should be very large and well shielded from back- 
ground radiations. The are to have either S > 1 km 2 , especially if sec- 
ondary muons (say from z/ M + iV — > /i~ + X) are detected, and/or gi- 
gantic mass, M > 1 x 10 9 ton, in the case of registration using nuclear- 
electromagnetic cascades (from hadron state X in the case or from 
v e + N — > X). These requirements may be met: 

• by search for Cerenkov radiation in deep underwater detectors with 
desirable volume under control V ~ 1 km 3 (mass M ~ 10 9 ton); 

• by registration of air nitrogen fluorescence, induced by Extended Air 
Showers; an effective control over M ~ 10 11 ton of atmosphere may 
be achieved either from the earth (see Fly's Eye || and HiRes 0) 
or from satellites (Airwatch || etc); 

• using such a 'neutrino' detector, as Pierre Auger installation^] |J. 
Note, that Pierre Auger telescope is designed for the study of CRs at 
E > lx 10 20 eV, high-energy neutrino physics being just a byproduct. 

1.4 Rate of events 

The rate of nucleon-electromagnetic cascades in a neutrino detector is pro- 
portional to convolution product of z/iV-cross-sections and z/-fluxes, and 
also to the number of scattering centers, i.e. mass of detector. A standard 
requirement is to have at least 10 events per year, but the only way to 

1 A hybrid detector in Argentina will consist of 1600 tanks, each filled with 12 m 3 
of water. Stations will be distributed in a grid with 1.5 km spacing, b) Four Fly's Eye 
type fluorescence detectors, controlling ~ 3000 km 2 of the site. 
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gain higher statistics is to increase the installation mass. But huge arrays 
are very expensive. However, the higher are i/iV-cross-section, the higher 
is the expected rate of events, the better results may be obtained. 



1.5 Background 

A severe problem for HENA is high CR background. One is to shield 
a detector either burying it deep underground (underwater, underice) or 
selecting just down-up going events, using the Earth as a shield. In any 
case, the atmospheric neutrino background (mostly from prompt z/'s) is 
inevitable. Nevertheless, many //-source models suggest that extragalactic 
diffuse z/-flux dominates at E > 10 15 -j- 10 17 eV. 



1.6 The aims of the paper 

We shall discuss several consequences for HENA of accelerated growth of 
neutrino-nucleon cross-sections at extremely high energies. We argue that 
such acceleration may be induced by non-perturbative hard pomeron with 
intercept ~ 1.4, which was proposed and successfully exploited in series 
of recent papers by A. Donnachie and P. V. Landshoff |DJ. We study 
its influence on the rate of cascades in a //-detector and on evolution of 
neutrino spectra during their pass through large column depths of matter. 



2 Cross-sections 

Neutrino-nucleon scatterings proceed via CC and NC interactions: 

vM+N^ F + X, (CC) (1) 

u&fi+N^ vi{vi)+X, (NC) (2) 

where I = e, /x, r. High energy behavior of z/iV-cross-sections is yet un- 
known. All attempts to estimate them at E > 1 TeV are mostly based on 
different extrapolations of nucleon Structure Functions (SFs) F£ N (x, Q 2 ), 
F» N (x,Q 2 ) to x < 1 x KT 5 (see e.g. Ref.s [0, [□], |lj§ 0). In papers 



[15|1 a new parameterization of these SFs has been derived. It was 
grounded on a successful description of F^(x, Q 2 ) by A. Donnachie and P. 
V. Landshoff ]17J , which claimed that small-x ep-scattering data of HERA 



1~6| may be successfully explained with the help of a simple combination 
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of several Regge theory inspired leading pole terms. The most important 
were so-called 'soft' (with intercept ~ 1.08) and 'hard 7 (with intercept 
~ 1.4, hence very similar to the perturbative BFKL one pomerons; 
first prevails at small Q 2 , while the latter dominates at large Q 2 . More- 
over, DL argued that perturbative QCD (pQCD) fails at small x < 1CT 5 
and even that its validity at x ~ 10~ 4 -j- 10~ 5 is a pure fluke. However, the 
DL approach is also just a model, neglecting, in particular, all non-leading 
poles and cuts in the complex angular momentum /-plane. This results in 
violation of unitarity and unlimited growth of cross-sections at E v — > oo. 

Using generalization of DL F^{x, Q 2 ) description to i/iV-scattering 
case, SFs F% N (x, Q 2 ) and F^ N (x, Q 2 ), presumably valid in the whole range 
of kinematic variables < x < 1 and < Q 2 < oo, have been derived in 
Ref. ||15|| . At x > 10~ 5 they were determined by pQCD parameterization 
CTEQ5 [TJ5|], while in the small- a; region they were driven by the analogous 
Regge theory inspired description. A special interpolation procedure have 
allowed to meet smoothly these different, both over x and Q 2 , descrip- 



tions of SFs at low and high x. In Ref. [|14| these new SFs were denoted 
by DL+CTEQ5 indicating that they had their origin in the interpolation 
between DL and pQCD descriptions. 

In parallel there were considered in Ref. |15| vN SFs, which had been 
derived by a simple extrapolation of pQCD SFs from x > 1 x 10 -5 to the 
small-x region. These were given by formula: 



pvN,Log+CTEQ5 ( _ ^ rf2\ 



[X < x r 



Q z 



^vN,CTEQb , 



X r 



,Q 2 



X 



(3) 



d\nF» N > CTEQ5 (x,Q 2 



dlnx 



X r 



1 x 10" 



(4) 



i — 2, 3. These SFs smoothly shot to the low- a; region from the CTEQ5 
defined high-x one [0, [14|, Starting values of functions and of their 
logarithmic derivatives over In a; had been taken at the x = x m i n bound- 
ary of CTEQ5. This parameterization was designated by Log+CTEQ5. 
Actually, these SFs are very close to analogous of Ref. JTTJ . 



Denoted as DL+CTEQ5, total CC + NC z/A^-cross-sections are shown 
in Fig. |T|. For comparison there are also plotted the corresponding cross- 
sections, which were obtained 



a) in the framework of simple Log + CTEQ5 extrapolation (||D; 
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b) with the help of CTEQ4 parameterization by Gandhi et al. Ref. 

denoted as GQRS-98 (CTEQ4); 

c) within a united BFKL/DGLAP approach by Kwiecinski, Martin and 



Stasto 13 , labelled as KMS. 



DL+CTEQ5 
Log+CTEQ5 
GQRS-98 (CTEQ4) 
KMS 




DL+CTEQ5 
Log+CTEQ5 
GQRS-98 (CTEQ4) 
KMS 
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Figure 1: Total (CC+NC) uN- 
cross-sections. Labelling of the 
curves is described in the text. 



Figure 2: r(E u ) are ratios between 
cross-sections shown in Fig. [I] and 
the Log+CTEQ5 cross-section. 



Due to hard pomeron, our model predicts the most rapid growth of 
cross-sections at high energies. The differences between cited calculations 
become especially clear in Fig. 0, where each cross-section is divided by 
the cross-section of Log+CTEQ5. 



3 Rates of cascades in z/-detectors 

In this section we shall check the manifestations of pomerons in the rates 
of nuclear-electromagnetic cascades in deep underwater /underice detec- 
tors. To simplify the problem, we shall deal with an incident z/-flux energy 
spectrum, approximated by a simple power-law decreasing function, 

F v {E u ^) = AxE-^ +1 \ (5) 

where 7 is the index of an integral neutrino spectrum, F v {> E,^^ and 
A is the normalization coefficient. Most source models commonly assume 
11 < 7 < 2.1. The real generation //-spectrum should be cut-off at high 

2 Note that for power-law spectra j x F u (> E,j) = E x F„(E,j). 
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energies of course, though in top-down models this may occur only at 
very large energy. Moreover, as it will be shown later, the high-energy 
part of neutrino spectrum suffers from attenuation on passage through 
large column depths of matter. Nevertheless, for the purposes of present 
consideration all these details may be temporarily neglected. 

So, in the case of power- law decreasing spectra (||), differential and 
integral rates of cascade production in a detector may be calculated with 
the help of so-called demensionless differential, 

1 

v it? \ I a ■ y da vN {E h /y ) y) 

Z h [Eh,V= / d W j , (6) 



and integral, 



1 

Y h {E h , j)=j f duu^ x Z h f — , 7 ) , (7) 



o 



hadron moments |Z0|, . Here Eh is the energy of a hadron-electromagnetic 
cascade, y = E h /E u , and <t = G 2 F m^ / /'n is the normalization cross-section 
(for rriyy = 81 GeV <Jq = 1.09 x 10~ 34 cm 2 ). These rates in the case of 
CC-scattering ([[]) are 

dN h (E h )/dt = Zf c (£ ft , 7 )AW^(^)> (8) 
dN h (> E h )/dt = Y^ c {E h)1 )N N a nF v {> E h ), (9) 

where Nn is the number of nucleons in a detector, Q is the effective solid 
angle, opened for the neutrino flux. To take into account the i>N- and NC- 
interaction cases, one should just substitute the appropriate differential 
cross-sections in the Eq. (|[) for the CC cross-section. 

In the case of electron (anti) neutrino CC-scattering (H) the role of 
differential hadron moment belongs to the corresponding normalized CC- 
cross-section, a^ N (E v ) / a ; since here E h = E u , the i/ e -fTux in (HH) is to 
be taken at the energy of incident neutrino. 

According to Eq.s (§,0), hadron moments are sensitive to the high- 
energy parts of i/-spectra. As a consequence, hard pomeron effects look 
even more pronounced in these observables. To illustrate a common trend 
of hadron moments, the differential and integral uN and vN hadron mo- 
ments are plotted in Fig. [3] and Fig. [|, respectively, for DL+CTEQ5 pa- 
rameterization and 7 = 1.1, both for CC and NC. 
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Figure 3: Differential CC and NC Figure 4: Integral CC and NC 
hadron moments for 7 = 1.1. hadron moments for 7 = 1.1. 



A detailed discussion of dependencies of hadron moments on 7 can be 
found in Ref. IE2. 
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Figure 5: Ratios between CC 
differential hadron moments of 
DL+CTEQ5 and Log+CTEQ5 pa- 
rameterizations for 7 = 1.1. 
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Figure 6: Analogous to Fig. |5] 
ratios of integral moments. The 
curve denoted v e N describes the 
case of E h = E u . 



It is informative to plot ratios between corresponding CC hadron mo- 
ments, evaluated within DL+CTEQ5 and Log+CTEQ5 parameterizations, 
respectively. They are shown in Fig.s HH for 7 = 1.1. These figures cor- 
roborate our suggestion, that hadron moments are more sensitive to the 
non-perturbative pomeron effects. The strongest effect in integral mo- 
ments takes place for v e {v e )N CC-scattering, when the whole energy of 
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the incident neutrino is transferred into cascade. 



4 Absorption and regeneration of neutrino 
fluxes in matter 

Neutrinos are known to be weak interacting particles. However, HE neu- 
trinos should be absorbed in large dense samples of matter, when a v ^{E) x 
n xl/rriN > 1; here n (in g/cm 3 ) is the matter density and I (in cm) is the 
range of neutrino pass. An appropriate variable for such consideration is 
column depth, r = n x I (measured in g/cm 2 ). 

Attenuation effects may be described by a shadow factor, which is a 
ratio of neutrino flux at column depth r and the primary flux (i.e. at 
r = 0): 

To a first approximation, shadow factor S(E, r) fa exp(—a cc+N (E)t /mw) 
It provides an estimate that neutrinos with E > 10 TeV (see Fig. |l|), pass- 
ing through the center of the Earth (column depth of matter in this case 
is r ~ 1.1 x 10 10 g/cm 2 ), are partially absorbed. 

On the one hand, the higher is the energy, the stronger is attenuation. 
On the other hand, Earth is the natural, but not the unique example of 
a large column depth of matter in the Universe. One can easily imagine 
different environments with much higher column depths, where neutrino 
flux attenuations would be stronger as well. Hence, attenuation effects 
should be taken into account in HENA. 

In this section we emphasize that hard pomeron enhanced i/iV-cross- 
sections bring to much higher attenuations of neutrino fluxes, than it was 
previously expected. 

The absorption of neutrinos in the Earth have been earlier discussed 
by many authors. However, only in paper |10| the regeneration of neutrino 



fluxes due to iVC-scattering (|2|) has been taken into account for the first 
time. To evaluate the effect, a 'thin target' model have been used in that 
paper. That model assumed that average attenuation length in the Earth 
is comparable with its diameter. And that was the case for the Earth as a 
target and the energy range considered in Ref. fTofl . 



Nevertheless, in order to treat correctly the higher energies and larger 
column depths, a more sophisticated consideration is necessary. Ten years 
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later it was proposed by A. Nicolaidis and A. Taramopoulos in Ref. I fZSfl . 
They have derived a trivial transport equation 



&F V (E V , r) = -a cc+NC {E u ) x F V (E V , r 
j_ r^jL. _ 

-y ay 



+ Jo ft >< d ° NC(E t- yly) >< -V),r) [ll > 



and solved it (unfortunately, under some wrong assumptions about differ- 
ential cross-sections). First term in the r.h.s. of the Eq. fllTD accounts for 
the absorption, while the second one describes regeneration via NC. 

Later these equations have been solved correctly by many authors 
with various parameterizations of high-energy z/iV-cross-sections (see e.g. 
Ref.s 0,11). 



An elegant method for solution of the transport equation (JTTJ) , along 
with a detailed discussion of different neutrino spectra evolution, has been 
given by V. A. Naumov and L. Perrone in paper |[24j| . In this paper we 



exploited their method, applying it to the infinite incident power-law de- 
creasing neutrino spectra (|5[). Calculations were performed for the two 
parameterizations of differential cross-sections, viz. for DL + CTEQb and 
Log + CTEQ5. 

Our results are demonstrated in Fig.s 0-|lO[ 
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Figure 7: Shadow factors for 
uN- and z/iV-scattering for sev- 
eral values of column depth, pa- 
rameterized by the factor k = 
t/t : t = 1.1 x 10 lc Wcm 2 cor- 

responds to the passage through 
the center of the Earth. 




2 4 6 8 10 12 



Log 10 E [GeV] 

Figure 8: Ratios of vN- 
scattering shadow factors 
calculated with and without re- 
gard for regeneration of ^-fluxes 
by neutral currents. Factors k 
are the same as in Fig. [7|. 
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In Fig. |7] shadow factors S(E, r) (|IUD, evaluated with the help of hard 
pomeron enhanced parameterization DL + CTEQ5, are plotted for neutri- 
nos and antineutrinos. All curves relate to only one value of initial neutrino 
spectrum index, 7 = 1.1. Evolution is illustrated by successively increasing 
values of column depth r, which are expressed in terms of the 'standard' 
column depth r = 1.1 x 10 10 g/cm 2 . 

In Fig. |g the importance of the regeneration term is emphasized. The 
ratios of shadow factors with regeneration and those accounting only for 
absorption are shown here for the same successively increasing column 
depths t as in Fig. [7|. 




1 2 3 4 5 6 7 8 9 10 11 12 



Log 10 E [GeV] 

Figure 9: Ratios of vN shadow fac- 
tors for incident z/-spectra with in- 
dices 7 = 2.1 and 7 = 1.1, respec- 
tively. Factors k are the same as in 
Fig-0. 




Figure 10: Ratios of vN shadow 
factors, calculated within DL + 
CTEQ5 and Log+CTEQ5 param- 
eterizations, respectively. Factors 
k are the same as in Fig. [7|. 



The dependencies of shadow factors on 7 may be understood from 
Fig. |9|. The ratios of these factors are plotted here for two indices of initial 
spectra, viz. 7 = 2.1 and 7 = 1.1. Evolution is shown in the same way as 
it was done in the previous figures. 

And, finally, the ratios of shadow factors calculated within DL+CTEQ5 
and Log + CTEQ5 parameterizations are shown in Fig. [10]. Notations for 
evolution remain the same as in the previous pictures. This plot persuades 
that at high energies the hard pomeron effects cannot be neglected even 
for rather shallow column depths of matter. Pomeron effects, being raised 
in shadow factors into exponent, bring to tremendous changes as compared 
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with previous estimations. These coefficients are evidently the most sensi- 
tive observables in the HE neutrino astrophysics to the hypothetical hard 
pomeron enhancement of vN structure functions. 

Conclusions 

In this paper we have studied different manifestations of hypothetical hard 
pomeron enhanced vN structure functions in the high-energy neutrino as- 
trophysics. We have compared i) cross-sections <j u n(E), ii) rates of hadron- 
electromagnetic cascades, parametrized through differential and integral 
hadron moments, Z h (E h ,^),Y h (E h ,^), and iii) shadow factors S(E v ,t), 
calculated within hard pomeron enhanced DL+CTEQ5 parameterization, 
with analogous values calculated within simple extrapolation of structure 
functions from high-x perturbative QCD description to the low-x region, 
Log + CTEQ5. 

We found that these effects, if our assumption are justified, may change 
substantially our notion of high energy neutrino astrophysics. The most 
pronounced evidence for their presence has been found in the rapid decrease 
of shadow factors with energy. It may be seen even with column depths 
available at the Earth. We hope the validity of our suggestions will be 
checked by future giant neutrino telescopes. 
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